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Sonoluminescence is explained in terms of quantum radiation by moving interfaces 
between media of different polarizability. In a stationary dielectric the zero-point 
fluctuations of the electromagnetic field excite virtual two-photon states which become 



Oh 



real under perturbation due to motion of the dielectric. The sonoluminescent bubble 
is modelled as an optically empty cavity in a homogeneous dielectric. The problem of 
the photon emission by a cavity of time-dependent radius is handled in a Hamiltonian 
formalism which is dealt with perturbatively up to first order in the velocity of the 
bubble surface over the speed of light. A parameter-dependence of the zero-order 
Hamiltonian in addition to the first-order perturbation calls for a new perturbative 
method combining standard perturbation theory with an adiabatic approximation. In 
this way the transition amplitude from the vacuum into a two-photon state is obtained, 
and expressions for the single-photon spectrum and the total energy radiated during 
one flash are given both in full and in the short-wavelengths approximation when 
the bubble is larger than the wavelengths of the emitted light. A model profile is 
assumed for the time-dependence of the bubble during the collapse, and in this model 
the radiated energy and the spectrum are calculated numerically and in the short- 
wavelengths limit also analytically. It is shown analytically that the spectral density 
has the same frequency-dependence as black-body radiation; this is purely an effect 
of correlated quantum fluctuations at zero temperature. The present theory clarifies 
a number of hitherto unsolved problems and suggests explanations for several more. 
Possible experiments that discriminate this from other theories of sonoluminescence 
are proposed. 
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1 INTRODUCTION 

1.1 State of the art 

Sonoluminescence is the phenomenon of light emission by sound-driven gas bubbles in fluids, 
ordinarily air bubbles in water. Sound makes bubbles collapse or expand, and a rapid flash 
of light is observed after each collapse. This phenomenon has been known for 60 years [T|, 
but came under systematic investigation only recently when experimentalists succeeded in 
trapping bubbles and maintaining sonoluminescence as a stable process over hours or even 
days 0,1. 

During stable sonoluminescence || H a bubble is trapped at the pressure antinode of a 
standing sound wave, which typically has a frequency of about 25 kHz. With an astonishing 
clocklike precision the bubble sends off one sharp flash of light per acoustic cycle. Less 
than 10 ps is commonly given as a conservative estimate of the pulse length. The observed 
jitter has been found to be extremely small and to show curious phase properties whose 
origin could so far not be identified Q. The spectral density of the light emitted drops 
with wavelength and resembles the tail of a black-body spectrum of several tens of thousand 
Kelvin [|. 

Whereas the dynamics of the bubble motion has been successfully explained and a theo- 
retical model by Lofstedt, Barber, and Putterman || based on rather involved hydrodynamic 
calculations reproduces the experimentally measured time-dependence of the bubble radius 
remarkably well, the process of the light emission has so far defied any theoretical eluci- 
dation. That is why the present paper focuses on the radiation process, making use of the 
knowledge about the hydrodynamics of the bubble motion as input. 

There have been several attempts of explaining the light seen in sonoluminescence. The 
apparent similarity of the spectrum to a thermal spectrum has led to the hypothesis that the 
light might come from a process of black-body radiation or bremsstrahlung [||, § . Along this 
line it has been argued that the gas in the collapsing bubble is compressed so strongly that a 
plasma is formed which then radiates. However, one can quickly convince oneself that neither 
black-body radiation nor bremsstrahlung can possibly account for the radiation observed in 
sonoluminescence. Either of them would lead to a continuous spectrum whose major part 
would lie below the absorption edge of water at 180nm and would therefore be absorbed by 
it. Estimating from the corresponding visible part of the spectrum the amount of energy 
that would be absorbed, one obtains such a large number that one would expect to see 
rather obvious macroscopic consequences of the absorption 0, as for instance dissociation 
of the water molecules, formation of radicals etc., which have not been observed. Moreover, 
black-body radiation is an equilibrium phenomenon and involves several atomic transitions; 
it could very unlikely explain pulse lengths of less than lOps. Neither is any explanation 
involving bremsstrahlung satisfactory, because it would entail the presence of free electrons 
and rather slow recombination radiation. 

Rather more convincing is Suslick's theory [TI| which explains the sonoluminescence spec- 
tra on the basis of pressure-broadened rotational and vibrational lines in diatomic emission 
spectra. For silicone-oil sonoluminescence one finds an excellent agreement of synthetic and 



observed spectra by considering emission from excited-state C2 [10|. For water, however, any 



attempts to model the spectrum on the basis of this theory have so far been unsuccessful, 



Sonoluminescence 



3 



although the well-known 310nm system of OH is thought to be largely responsible for the 
broad peak around this wavelength in the observed spectrum ||, [TT| . 



The most recent speculation on the nature of sonoluminescence radiation is a theory of 
collision- induced emission |12) , which, however, in its present version contains still too many 
indeterminate points and adjustable parameters to permit a judgement on its tenability. 

1.2 Quantum vacuum radiation as a candidate 

This article pursues a line of thought loosely inspired by Schwinger |T3j who suggested that 
sonoluminescence could be some kind of dynamic Casimir effect, which the present writer 
agrees with in so far as the light emission observed in sonoluminescence has, just like the 
Casimir effect, its origin in the interaction of the vacuum fluctuations of the quantized 
electromagnetic field with a dielectric medium. Sonoluminescence is, however, much more 
closely related to the Unruh effect. 



Let first the Casimir effect [14], |15| be recalled: two parallel conducting or dielectric 
plates in vacuum feel an attractive force which arises from the boundary conditions the 
plates impose on the vacuum electromagnetic field. In a more intuitive picture one can 
understand the Casimir effect in terms of van-der-Waals forces; the electromagnetic zero- 
point fluctuations induce local fluctuating dipoles in each of the plates and because of the 
spatial correlations of the fluctuations the interaction of these dipoles leads to a net attractive 
force. 

The Unruh effect [|T6| , [17| is a dynamic generalization of the Casimir effect and predicts 
radiation by non-inertially moving mirrors. This phenomenon is not exclusive to perfect 
mirrors, i.e. perfect conductors; quantum radiation by moving dielectrics has also been in- 
vestigated |18| |, and moreover some of the pathological points of the perfect-reflector theories 



can be circumvented in the more physical case of dielectrics. Again, the intuitive picture 
of the process is that the zero-point electromagnetic field excites fluctuating dipoles in the 
(perfect or imperfect) mirror and these dipoles are the source of radiation when the mirror 
moves non-uniformly. 

A more rigorous way of understanding why a moving mirror that interacts with the 
vacuum fluctuations of the quantized photon field emits radiation, is to start by considering 
a nominally stationary mirror, with the intention of eventually applying the fluctuation- 
dissipation theorem. The radiation pressure on the mirror is given by the vacuum expectation 
value of the force operator J 7 , which is obtained from the stress-energy-momentum tensor 
of the electromagnetic field subjected to appropriate boundary conditions on the surface of 
the mirror. The net force on a single stationary mirror in vacuum is of course zero by virtue 
of translation invariance, 

F = (0|J"|0) = . 

However, the mean-square deviation of this force does not vanish, since the force operator T 
does not commute with the Hamiltonian. In other words, the mirror is exposed to radiation- 
pressure fluctuations, whose mean-square deviation is given by 

AF 2 = (0|F 2 |0) - (0|F|0) 2 . 

Knowing that the force operator T is (just like the Hamiltonian) a functional that is quadratic 
in the field operators, i.e. quadratic in the photon annihilation and creation operators, one 
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can use the decomposition of the identity into projection operators onto a complete set of 
photon eigenstates, of which then only two-photon states survive, and rewrite fl~£fl : 



AF 2 = ^ J dk J dk' \(0\F\k,k'}\ 2 . (1.1) 



This means that virtual two-photon states are perpetually excited by the mirror in the 
vacuum, in accordance with the fluctuating radiation pressure. Yet the fluctuating forces 
on the left and right sides of the mirror are balanced against each other, so that no mean 
radiation pressure acts on the mirror. By virtue of Lorentz invariance, the same is true for 
a mirror that moves with constant velocity. 

However, when the mirror moves non-inertially, the radiation-pressure fluctuations on 
opposite sites of the mirror are out of balance and the mirror experiences a non-vanishing 
frictional force. The virtual two-photon states turn into real states, and the loss of momen- 
tum by the radiation of the photon pairs provides the physical mechanism for the friction felt 
by the mirror. The fluctuation-dissipation theorem puts this into formulae and interrelates 
the power spectrum of the fluctuations on the stationary mirror and the dissipative part of 
the response function that connects the force on the moving mirror to its velocity [^0 . 



It is a well-established fact that radiation by moving mirrors shows thermal properties 
although one is dealing with zero-temperature quantum field theory. The original statement 
of the Unruh effect [^, |I7|] is that a mirror moving with constant proper acceleration a 



in vacuum appears to be radiating particles as if it were a black body at a temperature 
Tunruh = ha / (2tt k^c) . The reason for this behaviour is that the photons are radiated in 
correlated pairs, in the language of quantum optics — they form a two-mode squeezed state, 
and the observation of the single-photon spectrum involves a tracing over the other photon 
of the pair which is well-known to entail thermal properties of the state [[HJ. Formally this 
connection is established by representing the two-mode state in a dual Hilbert space and 
making contact with the theory of thermofield dynamics 122 . 



As to an experimental verification of the Unruh effect, the record is empty. Understand- 
ably so, because the Unruh temperature is tiny for commonly achievable accelerations. The 
only viable suggestion for an experiment has come from Yablonovitch |23|] , who thought that 
the sudden ionization of a gas or a semiconductor crystal might produce an accelerating 
discontinuity in the refractive index fast enough to radiate a measurable amount of photons. 

From all of the above, quantum vacuum radiation seems to be a good candidate for 
explaining the radiation process in sonoluminescence. The surface of the bubble is the 
moving interface of discontinuous polarizability, i.e. the moving mirror. In the visible range 
water has a refractive index of 1.3, and the gas inside the bubble has a refractive index 
of practically 1 even when strongly compressed. Although the discontinuity of 0.3 in the 
refractive index is not huge, it is large enough to radiate an appreciable number of photons 
if the motion is sufficiently fast. In fact, the discontinuity in the refractive index will enter 
the final results for the radiated spectrum merely as a prefactor, and hence only its order of 
magnitude is important. 

Of much greater significance is the highly non-linear dynamics of the bubble motion. At 
the point when the bubble collapses and starts re-expanding, the velocity of the interface 
changes its sign. As known from experiments as well as from model calculations of the 
bubble dynamics || |S| this turn-around is extremely fast, which means that tremendous 
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accelerations and higher moments of the motion are involved. The present theory predicts 
a burst of photons as a consequence. 

As it will be shown in more detail in the course of this article, the theory of quantum 
vacuum radiation resolves several to date unexplained issues. 

The fact that the photons are radiated in correlated pairs leads to thermal properties of 
the one-photon spectrum, unrelated to the temperature in the bubble which is presumably 
far too small to cause any major effect. 

In accordance with a weakly frequency- dependent refractive index the radiation spectrum 
shows features at the resonance frequencies of the dipolar molecules in the medium. As water 
molecules are highly polarizable, this is to be expected a discernible effect; it explains the 
relation of the peak around 310nm in the spectrum to the well-known OH line. 

Barely any photons are created below the absorption edge of water, as the polarizability 
is far too small in this region. Therefore, few photons are absorbed and no macroscopically 
noticeable changes of the water are to be expected. 

The pulse length predicted by the theory of quantum vacuum radiation is of the order of 
the time it takes for the zero-point fluctuations to correlate around the bubble. With bubble 
sizes of around 1/im or less, the time light takes to cross the bubble is in the femtosecond 
range. Otherwise the time-scale is of course influenced by the dynamics of the motion of the 
bubble interface at and just after the collapse. 



1.3 Outline and overview 

The theory of quantum vacuum radiation by a gas bubble in water will be expounded in the 
following sections. Water will be understood as a non-absorbing dielectric describable by a 
constant refractive index. This is a good approximation in the spectral region of interest 
where water is only weakly dispersive. By virtue of adiabaticity the refractive index n can 
be replaced by n(u) in the end result for the radiated spectrum. The gas inside the bubble is 
optically so thin, even at the collapse of the bubble, that its refractive index will be assumed 
to be 1 throughout the calculation. 

The bubble will be considered as externally driven, i.e. the radius of the bubble as a 
prescribed function of time; the hydrodynamics of the bubble motion is not the concern of 
this paper. However, the back-reaction of the radiation process on the motion of the bubble 
will be specified. 

Hence the problem is reduced to a model of a spherical cavity of radius R(t) in a homo- 
geneous non- dispersive dielectric described by a constant refractive index n. The radiated 
spectral density will be obtained as a functional of R(t). 

The next section deals with the quantization of the photon field in the presence of a 
stationary spherical bubble in a homogeneous dielectric. In section III the Schrodinger equa- 
tion for the photon state-vector is written down, and the vacuum-to-two-photon transition- 
amplitude is calculated by a method of time-dependent perturbation theory that accommo- 
dates both an adiabatic time-dependence of the Hamiltonian and a perturbative addition to 
the zero-order Hamiltonian |24| . Section IV states and examines the results for the radiated 



energy and the spectral density. The appearance of a thermal-like spectrum is demonstrated 
and numerical results are presented. Finally, section V gives a summary and a critical re- 
flection on the theory of quantum vacuum radiation by sonoluminescent bubbles. Strengths 
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and weaknesses of the theory are scrutinized, and open questions are voiced. 

Readers not interested in the technicalities of the theory are encouraged to look over 
section IV for an aggregate of the essential results and to read section V for a guided summary 
and interpretation of all results. 

Several appendices contribute technical details necessary for the clarity of the presenta- 
tion. Appendix A calculates the Hamiltonian for a uniformly moving dielectric in preparation 
for section III. Appendix B gives the mode-expansion for the Helmholtz equation in spherical 
coordinates, which is essential throughout the paper. The force on a stationary dielectric is 
determined in appendix C. 

CGS units are used everywhere in the paper; h and c are set equal to 1 unless explicitly 



indicated. All special functions are defined as in refs. [25, 26 



2 QUANTIZATION OF THE PHOTON FIELD 

The Hamiltonian for the electromagnetic field in the presence of a medium with dielectric 
function e(r) reads 

H, = \ fdh (- + B 2 ) . (2.1) 



2 J \ e 
A bubble of radius R is described by 

e(r;R) = 1 + (n 2 - 1) 0(r - R) , (2.2) 

where 9 is the Heaviside step function. This is to say that the dielectric constant equals 1 in 
the interior of the bubble and n 2 in the surrounding medium. The Maxwell equations imply 
continuity conditions for the fields across the boundary; these are: 



, D ll 

and — continuous 

e 

B continuous , 



(2.3) 



or in spherical coordinates 



T-joutside j-joutside 

"pvinside "pjoutside "pvinside 9 -pvinside 4> 

YJ — U , Uq — , YJ± — , 

rr rr (2.4) 

B inside td outside "D inside "D outside td inside "D outside 

r — D r } °6 — D 9 5 D 6 ~ D d> 



The Hamiltonian ( |2.1| ) depends parametrically on the bubble radius R via the dielectric 
function e(r; R). Although a problem with a varying bubble size is aspired to be solved, for 
the purpose of quantizing the photon field the radius of the bubble will be kept constant. In 
order to quantize the system for a time-dependent radius R(t) one would need to know the 
eigenfunctions of the time-dependent Hamiltonian ( |A6|) ; but knowing them would amount to 
the exact solution of the whole problem which is of course unachievable. As the calculation 
to follow in the next section will employ perturbation theory to first order in the velocity of 
the bubble surface (3 = R(t) over the speed of light in vacuum, a quantization for constant 
R is fully sufficient. 
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The field is quantized by ascribing operator nature to the field variables and imposing 
canonical commutation relations for the vector field A and its conjugate momentum II = 
— D. These are most easily implemented by expanding the field operators in terms of photon 
annihilation and creation operators, a k and a k respectively for a mode of momentum k and 
polarization s, and demanding that the latter fulfil the standard commutation relations 

" s s' 

a k , a k , 
a k , a k , 



= 6(k - k') 5 S 
. 



(2.5) 



At the same time the normal-mode expansions should be chosen such as to diagonalize 



the Hamiltonian (2.1) to the Hamiltonian of a photon field 



(2.6) 



All this is achieved by the following expansion of the electric displacement D and the 
magnetic field B: 



D TE 


= sfd*k *> [ 


B TE 


= v ^/d 3 k^ 


D TM 


= efd*k % 


Btm 


Juj 



a™ A™ + H C 
al M A™ - H C 



a™ A™ + H C 



(2.7) 



The mode functions An i2 ) are two linearly independent solutions of the Helmholtz equation. 
They satisfy the Coulomb gauge condition V ■ A( 1)2 ) = 0. The fields have been decomposed 
into their two transverse polarizations, chosen in spherical coordinates as the transverse 
electric (TE), for which the radial component of the displacement D vanishes, and the 
transverse magnetic (TM), for which the radial component of the magnetic field B vanishes. 
The mode functions A(i i2 ) and their properties are spelled out in appendix B. 

As mentioned above the quantization procedure is performed at an arbitrary but constant 
bubble radius R. This implies that although the Hilbert space of the quantized system 
stays always the same, the set of base vectors spanning it changes with R. A unitary 
transformation from the base at radius R to the one at radius R 1 exists in principle, but is 
hard if not impossible to find explicitly. So the bubble radius R serves as a parameter in 
traversing a continuous sequence of bases, and in a strict notation the photon annihilation 
and creation operators and the photon eigenstates should be supplemented by a label R. 
The vacuum or ground state of the field is defined by 

a s k (R) \ 0; R) = , (0; R | 0; R) = 1 ; 

single-photon states are written as 

| k s ; R) = a s k \R) | 0; R) , (k s ; R | k' s ,; R) = 8(k - k') S ss > ; 
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two-photon states are denoted by 

| k s , k' s ,; R) = a k \R) ap(R) | 0; R) , 

(k a , k' s ,-R\ l p , l' pl ; R) = 8(k - I) 5 sp 5(k' - I') 5 s , p , + 8(k - I') 5 sp , 5(k' - I) 5 s , p ; 
and so on for all higher photon number states. 

3 TWO-PHOTON EMISSION IN FIRST-ORDER PER- 
TURBATION THEORY 

The evolution of the state vector of the photon field is governed by the Schrodinger 
equation 

i±\i/>) = [H Q {R) + AH(R,P)]\1>) (3-1) 
where, according to the Hamiltonian (|A6|) derived in appendix A, 



H = \J d 3 r ^ + B 2 ) , (3.2a) 
AH = (3 f d 3 r £ -^- (D A B) r , (3.2b) 



and j3 = R is the velocity of the bubble surface. 

For the present purposes antisymmetrization of the operator product in AH can be 
dispensed with. What will be extracted from the mode expansion of (D A B) r are products 
of two photon creation operators a\a) k , which induce two-photon transitions from the vacuum. 
Since, however, creation operators commute mutually, operator ordering is inessential. 

To describe the sonoluminescence process by the Hamiltonian (3.2) means to ignore 
variations of the refractive index due to the periodic compression of the water in the vicinity 
of the bubble. This is a crude but innocuous approximation as long as the energies of the 
phonons excited in the water stay below those of the emitted photons. 

Initially the photon field is in its vacuum state while the bubble is at rest and has some 
radius R. As discussed at the end of the preceding section, the photon eigenstates depend 
parametrically on the radius R of the bubble. Hence the initial condition for the state vector 
\ip) reads 

|V(*o)> = |0;i*(to)>. (3-3) 



The integration of the Schrodinger equation ( |3.1|) poses a non-trivial problem since stan- 
dard methods of perturbation theory cannot be applied. The Hamiltonian AH ( |3.2b| ) cannot 
be treated as an ordinary perturbation because AH as well as H depend on the parameter 
R. The established way of dealing with slowly parameter-dependent Hamiltonians is the 
adiabatic approximation |27j]. However, the standard adiabatic approximation requires the 



knowledge of the complete set of eigenfunctions of the Hamiltonian for any allowed value 
of the parameter. In the present case only the eigenfunctions of part of the Hamiltonian, 
namely those of Hq, are known. Hence what is required is a judicious combination of standard 
perturbation theory and the standard adiabatic approximation; needed is a theory that is 
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capable of dealing both with a perturbative interaction Hamiltonian and with a Hamiltonian 
depending on a slowly varying parameter. 

Following the adiabatic theory by Pauli ||27|| , one starts with the eigenvalue equation for 
the unperturbed Hamiltonian Hq, solved for all possible values of the parameter R 

H (R) \n(R)) = E n (R) \n(R)) , (3.4) 

where E n (R) is the nth eigenvalue and \n(R)) the corresponding eigenvector. Here n is just 
a label; the eigenvalue spectrum need not be discrete. In general the levels can be multiply 
degenerate, so that |^(-R)) in fact stands for a whole subspace of orthonormal eigenvectors 
to the same eigenvalue E n (R). Where degeneracy matters it will be explicitly indicated by 
states \n'(R)) also belonging to E n (R). 

Differentiating ( |3.4j ) with respect to R and calculating the overlap with a state (m(R)\ 
one obtains 

/ I ^ I \ I 771 / I 3 I \ I I \ | T7T / I d I \ 

(ml— — \n) + E m (m\— — \n) = — — (m \n) + E n (m\— — \n) . 

Thus, provided no level-crossing occurs, ie if for m ^ n is E m (R) — E n (R) is different from 
zero for all possible R, as it will be the case in the present application, one has 

H-^» = ;g \ E H^rH for m^n. (3.5) 

Seeking a solution of the Schrodinger equation (|3.1|), one expands the wave- vector \ip(t)) 
into the eigenvectors of the instantaneous H (R(t)) 

\m) = t \n(R(t)))HR(t))\m) ■ (3-6) 
•fa 

Then the Schrodinger equation ( |3.1| ) becomes 

t i (4 |n> )^ <n ' V;)+i|n> (^ <n|V,> ) ^t 1 [^i n )H^) + Ai?|n}(n 
which by taking the scalar product with an eigenstate (m\ is turned into 

i{3(n\^)(m\-^\n) + iJ^(ra|V>) = E m (m\ip) + ^ (m\AH\n)(n\ip) . 

From here, application of the relation ( |3.5| ) yields 

d ^-A d 

i— (m\ip) -E m {m\ip) +i(3 ^ 

m'{E m ,=E m ) 



n(n^m) £ ' n £/ ™ ' ~ 



where the sum over m' takes states degenerate with m into account. Since eventually the 
transition probability Km^)! 2 and not the transition amplitude (m\ip) will be of physical 
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interest, one can gauge away the second term on the left-hand side of the above equation by 
defining 

(3.7) 



(m 



one finds 



den, 
dt 



Cm exp 



d 



-i / dr E m (r) 

J to 



t c m'(m\-^-\m') 

m'(E^=E m ) OU 



" t. 



(ml \n) exp 



i [ dr (E m - E ri 



E n — E m dR 
-i^ c n (m\AH\n) exp i^ dr (E m - E n ) 



(3.8) 



which is the key formula for the present approximation method. If it were not for the term 
containing AH, this expression would lead to merely the standard adiabatic approximation 
(cf. for instance ref. fl28|). It should be noted that both terms on the right-hand side of 
( |3.8| ) are of the same order, namely (3 1 ; the matrix element of OHq/OR is multiplied by j3, 
and AH is itself of order (3. 

The initial condition ( |3.3| ) for the wave function translates into the following initial 
conditions for the coefficients c m {t) defined by ( |3.7| ) 



c {t c 



c m ^o(^o) 







(3.9) 



As soon as the bubble starts moving, ie the interface velocity (3(t > t ) becomes different 
from zero, the rate of change of the c m is non-zero, as described by eq ( |3.8|) . For times t > to 
one has co(t) ~ 1 and c m ^o(t) = O (0) or higher. Hence, working only to first order in (3, 
one has to retain only the vacuum state in the summation over n on the right-hand side of 
eq (pp. 

Handling the time-dependence of the energy eigenvalues E n requires special care. For a 
cavity, like the bubble in the present problem, one has two limiting cases. The first is that 
the cavity walls are very poor reflectors; then it is convenient to label the cavity modes by 
wavenumber or energy as these are adiabatically conserved. However, if the cavity has a very 
high Q value, i.e. is close to perfectly reflecting, the number of nodes of the eigenfunction 
will be the adiabatically conserved quantity and not the wavenumber. The probability of 
reflection from an interface of a medium of refractive index n with the vacuum is given by 
(n — l) 2 / (n + l) 2 . For an air-water interface in the visible spectrum where n = 1.3 this leads 
to a reflection probability of less than 2%. Hence the bubble is a poor-quality cavity and 
the energy eigenvalues are adiabatically conserved. This justifies the use of wavenumbers for 
the labelling of eigenstates. 

Since H , and hence dH Q /dR, and AH are quadratic in the fields, the only transitions 
they can induce from the initial vacuum state lead to two-photon states, whence to first 
order in (3 p{| the system of differential equations (|3.8|) for the coefficients c m reduces to 



dc, 



kk' 



dt 



dH 0l 



uj + U) 



-(k,k';R\-^-\0;R)e 



j(u>+L>')(t-t ) 



-i(0; R\AH\0; R) , (3.10a) 
i(k, k'; R\AH\0; R) e i{uj+UJ ' ){t - to) . (3.10b) 
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The first equation is of no special interest, but the second equation will readily provide the 
solution of the problem posed — the perturbative description of the photon creation by the 
moving bubble surface. According to ( |3.7| ) the transition amplitude from the vacuum into a 
two-photon state is given by 

(k,k' ] R\^) = c kk! (t)e-^ + ^ t - t0 \ (3.11) 
and the initial vacuum evolves into the state 

l^> = |0; R) + n / dk dk' c w (t) e" 1 ^ >('-*>> \k, k'; R) . (3.12) 

Z J —oo J— oo 

The factor 1/2 takes care of the identical photon states \k,k';R) = \k',k;R), not to be 
double counted. 

It remains to evaluate the two matrix elements in eq ( |3.10b| ). With H as in ( |3.2a| ) one 
finds for the first 

(k,k';R\^-\0;R) = l - (l - R 2 j> dQ (k, k'; R\ D 2 r + n 2 E 2 \0;R) , (3.13) 

in the obvious notation En = {Eq^E^). Deriving this, one should bear in mind that H 
depends on R both through e(R) and through the discontinuity of Dn = (D e , D^). 
With AH as in ( p. 2b ) the second matrix element in ( 3. 10b| ) reads 

(k,k';R\ AfflO; R) = (l - -A I d 3 r (k, k'\ R\ (D A B) r |0; R) . (3.14) 

V n 2 ) Jr>R 

As photon states are eigenstates of Hq, one can write 

i (u + uj') (k, k'; R\ (D A B) r |0; R) = (k, k'; R\ i [H , (D A B) r ]|0; R) 

= (k, k'\ i?| ^ (D A B) r |0; R) . (3.15) 

to first order in (3. In further manipulating this expression, one can make use of the classical 
energy-momentum conservation law in a bulk dielectric 

^(DAB),, + V,T^ = 0, (3.16) 
where the stress-tensor in the medium is given by 

T ij = - B iBj + ^ + B 2 ) . (3.17) 

Since the photon field is a non-self-interacting field in Minkowski space, there is no doubt 
that this conservation law is valid also quantally. Using this and the above relation ( 3.15 ), 
one can rewrite ( 3.14 ) 

(k, k'\ R\ AH\0; R) = i^— (l - -V) / ( k i k ', R \ V/T^'lO; R) . 

UJ + UJ' \ n 2 / Jr>R 
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Applying Gauss' theorem leads to 



(k, k'- R\ A#|0; R) = -i^— (l--^) I dfi (k, k'- R\ T rr |0; R) 

to + to' \ n 2 J J 



( 1 - \) i dQ (k, k'- R\ I (n 2 E 2 -^ + B 2 B 2 r ) |0; R) . (3.18) 



12 



Inserting ( |3.13| ) and ( |3.18| ) into the expression ( |3.10b|) one finds for the time-derivative 
of the transition amplitude between the vacuum and a two-photon state \k, k') 



j^fdn (k, k'\ R\ (l + — ) D 2 r - Bj + B 2 r |0; i2) . 



<9t V n 2 J 2 to + to 

Comparison with the force operator T T given by eq. (|C5|) uncovers the relation of the photon 
creation to the pressure on the bubble; dc kk i / dt can be re-expressed as 

^ = ^— e 1 ^)^) (jfe, fc'; R\ T r 10; i?> . (3.19) 

at + a;' 

This is a truly remarkable result as it exposes the fluctuations of the radiation pressure as the 
origin of the photon-pair creation. As described by eq. flP ) for the mean-square deviation of 
the force, the fluctuations on a stationary mirror are tied in with the excitation of virtual two- 
photon states. The (non-uniform) motion of a mirror or a dielectric interface makes these 
virtual states become real, which has been explained in the two paragraphs following eq. 
( |1 . 1|) and is shown manifestly by eq. ( |3.19| ). The fluctuation-dissipation theorem underlies 



this connection; it however predicts only the dissipative force acting on the moving interface 
and not the photon-creation amplitude. Hence, the fluctuation-dissipation theorem cannot 
supersede the above derivation, but is nevertheless a useful check on it pi 



The integration of eq. ( |3.19 ) is complicated by the fact that the force matrix-element 
depends parametrically on R and therefore on time. Formally, the transition amplitude reads 

c kk ,(t) = l — f dr P(r) e^ + ^ T -^ (k,k';R{r)\ J>|0; i2(r)) . (3.20) 

00 + uo' Jto 

The matrix elements of the force operator, which is given by ( |U5|) , are calculated by ex- 
panding the D and B fields into normal modes as detailed in section 2 and appendix B. An 
unsophisticated calculation yields 

(k TE ,k' TE ;R\T r \0;R(r)) = (l - 1) £ Sj E 1 (k) Sj E ~\k') ^ 

x {£(£ + 1) Je (kR) Je (k'R) + [kR UkR)]' [k'R j £ (k'R)]'} 

x y™(k)Y7"*(k') , (3.21a) 

(k TM ,k' TM ;R\T r \0;R(r)) = (l - 1) V^7 £ $™ _1 (*) ^ ^ 

x [kk'R 2 - £(£ + 1)} 3e{kR) je{k'R) Y e m {k) Y e m *{k') , (3.21b) 
(k TM ,k' TM ;R\ F r \0;R(r)) = 0, (3.21c) 
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where the normalization constants iSj E,TM ~ 1 (k) are as given in eqs. ( |B8| ) and (|B9|). 

The probability of creating a photon pair in the mode \k, k 1 ) from the initial vacuum state 
is given by the modulus square of the amplitude Ckk', eq. fl3.20| ), featuring the above matrix 
elements of the force operator. The state of the photon field is specified by eq. (|3.12|) . Hence 
Ckk' carries all information one needs to determine the expectation values of all interesting 
observables, especially the total radiated energy and the spectral density, which are the 
subject of the following section. 



4 RADIATED ENERGY AND SPECTRAL DENSITY 

As I Ckk'(t)\ 2 1S the probability of creating a photon pair in the mode \k, k'), the total energy 
of the photons radiated by the bubble interface during one acoustic cycle reads 



i r T 

l JO J-oo 



oo roo 

d 3 k ' j31 -' 



d 3 k' (uJ + Uj') |Cfcfc/(t)| 5 



(4.1) 



where T is the period of the sound field. Inserting Ckk'(t) from ( 3.2(J| ) and (3.21) one obtains 
for W 



[n 



00 f°° uu} 

duj / da; , 

o Jo uj + uj' Jo 



T dr [ T dr' P(t)0(t') e^+^'XT-T') 
Jo 



x^(k,k',R(r),R(T')), (4.2) 



with the auxiliary function S defined as 

Z(k,k\R(T),R(T'))=J2W+V 



x {SJ E -\kR(r)) Sj E -\k'R(r)) Sj E * -\kR(r')) Sj E * ~\k' R(r')) 

x [£(£ + 1) j e (kR(r)) j e (k'R(r)) + \kR{r) je(kR(r))] f [k'R{r) j e (k'R{r))}' 
x ''£(£ + 1) h{kR{r')) h{k'R{r')) + [kR(r f ) ji(kR(r')))' [k'R(r') j e (k> 'R(r'))}' 

i 



+ S™ -\kR{r)) S™ -\k'R{r)) S™* -\kR{r')) S™* -\k> ' R{r')) 



x 



kk'RHr) -£(£+!) 



kk'R 2 '- 1 



T 



v. 



X j e (kR(T))j e (k'R(T))j e (kR(T'))j e (k'R(T'))} , 

and the normalization factors Sj E,TM ~ 1 (k) as given by eqs. 



(4.3) 



and (|BSp. 

A spectrometer measures the single-photon spectrum, which by symmetry is isotropic for 
a spherical bubble. The quantity of interest is therefore the angle-integrated spectral density 
radiated during one acoustic cycle, 



V(uj) 



UJ 



dt <i> dVt-L 



„L k I d s k' \ckk'{t)[ 

J J — oo 



(4.4) 



which becomes 

(n 2 - l) 2 



V(uj) = v " - ? uj 2 duj' , 

4:Tr z n z Jo {uj + uj y Jo 



UJ 



T dr [ T dr' P(t)(3(t') S"+"')^-t>) 
Jo 



x$s(k,h',R(T),R(T')) . (4.5) 
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The main difficulty in calculating the radiated energy W and the spectrum V(u) is the 
evaluation of the auxiliary function 3. To find an analytical approximation for 3 is a rather 
laborious task. 

First, note that both W and V(u), eqs. ( |4.2| ) and (|4.5|), contain a prefactor {n 2 — l) 2 
multiplying 3, so that they vanish in the limit n — > 1, as they should; in this limit there is 
no dielectric interface to produce radiation. This justifies an expansion around n = 1 in 3, 
which to first order does nothing but reduce the normalization factors Sj E '™ _1 in (|4.3|) to 
1, by virtue of eq. ( |B10| ); otherwise 9 is independent of n. 

In the present application the arguments of the Bessel functions in 3 are generally greater 
than 1, partly appreciably much greater, so that expansions for Bessel functions of small 
arguments are of no use here. As the summation over the index i runs up to infinity, 
the sum will be dominated by terms for which argument and index of each of the Bessel 
functions are comparable in magnitude. Hence Debye's uniform asymptotic expansion has 
to be employed (cf. [^] formulae 9.3.3 and 9.3.7). So, for instance, one obtains in the regime 
x>(£+l/2) 

jt(x) — ► - - : cos(z/tan/5 -ufi-j) 

v vsecptan/5 4 

and 

r ■ / \i/ /tan/3 7T. 

[x je{x)\ — ► -y— g sm(vtanp- vp - -) , 

where the abbreviations v = (l + 1/2) and x = i^sec (3 have been introduced. 

With the help of the above asymptotic approximations and by turning the summation 
over I into an integration one can derive that 3 behaves approximately like kk' 'R(t)R(t') 
in the short-wavelengths regime, i.e. when the photon wavelengths are shorter than the 
minimum bubble radius. Numerical investigation of the behaviour of 3 confirms this and 
yields 

3 ~ 1.16 kk' R(t)R(t') . (4.6) 

Employing this approximation one can integrate the expression ( |4.2| ) for W and obtains after 
a short calculation: 

(n 2 - I) 2 1 r T , d 5 R 2 (r) 



W = 1-16 ^ I dr R(r W r) . (4.7) 

One of the interesting consequences of this result is that the dissipative force acting on 
the moving dielectric interface can be seen to behave like R 2 (3^\t) (+ terms with lower 
derivatives of /3). This dependence tallies with results of calculations for frictional forces on 
moving perfect mirrors (see esp. [0]); the dissipative part of the radiation pressure on a 
moving dielectric or mirror is proportional to the fourth derivative of the velocity. 

The expression ( |4.7[ ) indicates also that any discontinuity in (3^ (t) or lower derivatives of 
(3 unavoidably leads to a divergence in W (and also in the spectral density V{uS)). Especially, 
one is not permitted to assume a step-function profile for R{t) during the collapse of the 
bubble, since this would give the physically meaningless result of infinite photon production, 
which is not salvageable by a cut-off or any other artificial regularization [13| . 



What produces the massive burst of photons from a collapsing sonoluminescent bubble 
is the turn-around of the velocity at the minimum radius of the bubble. There the velocity 
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rapidly changes sign, from collapse to re-expansion of the bubble. This means that the 
acceleration is peaked at this moment and so are higher derivatives of the velocity. 

In order to estimate the total energy radiated during one acoustic cycle, one can use the 
approximation (|4.6|) to rewrite the expression (|4.2|) as 



W= 1.16 



1; 



960 n 2 n 2 Jo 



dR 2 (r) 



o or 



(4- 



where Q is the sum of the photon frequencies in a pair. As W is a functional of the time- 
dependent radius R(t), one has to model R(t) appropriately in order to be able to obtain a 
number for W. At the collapse of the bubble the function R(t) has a sharp dip; hence it is 
reasonable to adopt the model profile 



R*(t) = Rq — (^Rq — R 



(t/7) 2 + 1 



(4.9) 



for the time-dependence of the bubble radius. Figure 1 illustrates R 2 {t) for various values 
of the parameter 7 which describes the timescale of the collapse and re-expansion process. 
The shorter 7 the faster is the turn-around of the velocity at minimum radius and the more 
violent is the collapse. In the figure 7 is half of the width of the dip halfway between R^ and 



i? min . In this simple model the total radiated energy (|4.8|) reads in SI units 



W = 1.16 



3(n 2 



h 



512 n 2 



c 4 7 5 



(2 2 \ * 

-^0 — -^min J 



(4.10) 



Experimental data on sonoluminescent bubbles @ suggest that Ro ~ 
0.5/iin are sensible values to assume. With n ~ 1.3 one obtains 



10/iin and i? E 



W = 1.8 ■ 10" 13 J for 7 ~ lfs. 



(4.11) 



which corresponds roughly to the experimentally observed amount of energy per burst. Cal- 
culating the radiated spectral density in the same model gives 



V{u) = 1.16 



in 



64 n 2 



c 4 7 v 



R min u e 



(4.12) 



This is one of the most important end-results of this calculation, as it exhibits the same us 



dependence as black-body radiation. Equating the exponent in (|4.12j ) to hcu/kT one derives 
that a turn-around time 7 of lfs corresponds to a temperature of around 4000K. This is 
however just a very crude estimate, as a lot of simplifications and approximations have been 
made in proceeding from ( f4.2p to ( ^-12j ); in general, the functional dependence of V on uj 
will not be as simple, although its overall behaviour is as characterized by eq. (|4.12 ). 

The major problem in the above derivation is that the approximation (|4.6j ) is good only 
at photon wavelengths that are smaller than the bubble radius R, i.e. when all products 
kR and k'R are greater than 1. Once kR reaches down to the order of 1 or even below, one 
has to expect resonances of the photon wavelengths with the bubble radius. An exploration 
of any such resonance effects requires taking into account the full kR dependence of the 
auxiliary integral 3 in eq. ( f4.3| ), which is fairly difficult even numerically. 
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The results of a computer simulation of a model like the one specified by ( |4.9| ) with 
Ro = 45/iin, -R m i n = 3/im, and 7 = lOfs are shown in figures 2 and 3. Comparing these 
with the values of V predicted by eq. ( |4.12| ), one sees that the numerical results show an 
enhancement of about a factor 1000 in V relative to the analytical approximation. This is due 
to resonant behaviour in (|4.9Q , since -R m i n is not any more appreciably much larger than the 
wavelength of the observed light. Numerical studies in the regime kR m i n ^5 1 are hindered 
by substantial expense in computation time; work by the present author is in progress [31fl. 
One can expect to see an even greater enhancement in V over the predictions of eq. (|4.12| ) 
for what are believed to be realistic values of -R m i n , i.e. values around 0.5/im. Thus, one can 
presumably substantially relax the requirement on how small the turn-around time 7 between 
collapse and re-expansion of the bubble has to be in order for the present theory to yield 
the experimentally observed number of photons. The crude model that led to the estimate 
( [4.1 1| ) would demand 7 to be as short as lfs, if it were to account for the experimental data 
for the sonoluminescence of an air bubble in water. The numerical calculation resulting in 
the figures 2 and 3 suggests that 7 ~ lOfs is fully sufficient; calculations for -R m i n smaller 
than 3/im will most likely require merely 7 ~ lOOfs. . . lOps, which is closer to what one 
expects this timescale to be on physical grounds. 

The features seen in the spectrum in figure 2 seem to be due to resonances between 
photon wavelengths and bubble size, but again, extensive numerical studies are needed to 
explore them |JT|. Comparison of the figures 2 and 3 makes it very clear that plotting 
data over photon wavelength rather than frequency tends to conceal such features; it might 
therefore be beneficial to plot experimental data over photon frequency as well. 

Given the time-dependence of the bubble radius R(t), one can in principle, technical 
difficulties aside, predict the sonoluminescence spectrum radiated by the bubble from eqs. 
( |1.5|) and ( f4.3|) . The next section summarizes the questions answered by the theory of 
quantum vacuum radiation and spells out some of the as yet unanswered ones. 



5 SUMMARY AND CRITICAL REFLECTION 

5.1 Successfully resolved issues 

5.1.1 Results for the spectral density and restrictions on the turn-around 
time 7 

The results of the previous section seem in concord with the experimentally observed facts 
in sonoluminescence. Knowing the time-dependence of the bubble radius R(t) one can 
evaluate the radiated spectral density from eqs. ( |4.5|) and ( |4.3| ). The analytical estimate 
( [4.1 1| ) for the total energy radiated during one acoustic cycle and the results of numerical 
calculations presented in figs. 2 and 3 show good qualitative and quantitative agreement 
with the experimental data ||. The function R(t) has been modelled by a dip, as written 
down in eq. ( |4.9| ) and made visual in fig. 1. The width of this dip is characterized by 
the important parameter 7, whose physical significance is that it is a measure of the time- 
scale of the turn-around of the velocity between the collapse and the re-expansion of the 
bubble. To model the experimentally observed data on the basis of the numerical calculations 
worked through so far, this turn-around has to happen in about lOfs. Further numerical 
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analysis of eqs. ( |4.5| ) and (O) in the regime where the radiated photon wavelengths get in 
resonance with the bubble size can be anticipated to relax this requirement to a turn-around 
time of roughly lOOfs. . .10ps, which seems quite realistic. Unfortunately, it is notoriously 
difficult to determine 7 experimentally, as it is arduous to measure the bubble radius close 
to the collapse; the laser light that is used to determine the bubble radius by means of 
fitting the scattering data to the Mie theory, is then most likely reflected by the shock- 
wave front propagating through the water rather then by the actual bubble surface 0. In 
addition, one might also have to take into account that the bubble shape deviates, perhaps 
even substantially, from spherical. Although this does not affect the present theory of the 
radiation mechanism to any discernible extent, since the spectrum of the vacuum fluctuations 
is known to be affected by the shape only to higher order [p3|, it will noticeably alter the light- 
scattering properties so that the Mie scattering theory is no longer applicable. Nevertheless, 
the currently available experimental data seem to suggest that 7 lies somewhere in the 
interval lOOfs. . .lOps0. 



5.1.2 Thermal properties of the spectral density 

The similarity of the observed photon spectrum to a black-body spectrum has its origin 
in the fact that the photon radiation emerges in coherent pairs. To obtain the single- 
photon spectrum, which is the one measured in spectral analyses, one has to trace over one 
photon in the pair, as done by integrating over k' in eq. ( [4.41 ). This process of tracing 
is known to engender thermal properties of the single-photon spectrum, even though the 
original two-photon state was a pure state and one has dealt with zero-temperature quantum 
field theory throughout |^TJ. In other words, it is the particular correlations within the 
radiated photon pairs that engender the thermal-like properties of the spectrum, and thermal 
processes are completely independent of this and in the case of sonoluminescence presumably 
of no significance whatsoever. 



5.1.3 Features in the spectrum 

Features in the otherwise smooth experimental spectra can be explained by a combination 
of three things: 

(i) Resonance effects in the radiation mechanism occur when the size of the bubble R(t) 
and the photon wavelengths A are of the same order of magnitude. The features seen in the 
numerically calculated spectral density in figure 2 are due to higher-order resonances of this 
kind where the R is an integer multiple of A; for direct resonances one expects much stronger 
effects. General predictions for such features can hardly be made as their detailed qualities 
are determined by the time-dependence of R, which leads over to the next point. 

(ii) Choosing gases other than air for the bubble contents leads to a modified dynamics of 
the bubble surface, as gas solubilities in water vary; but even a slightly different R(t) around 
the collapse will change the structure of the resonance effects between A and R, which is why 
one expects a strong dependence of the sonoluminescence spectra on the gas that saturates 



the water, as indeed observed in experiments [|TT], |32| . 

(Hi) The experimentally seen features in the spectrum might just as well be caused by 
dispersion, i.e. the dependence of the refractive index n on the photon frequency u. As long 
as n depends only weakly on u>, one can replace n by n(u) in eqs. Q4.5| ) and ( |4.12j ) by virtue of 
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adiabaticity. This explains that the experimental data show features at frequencies close to 
the vibration-rotation excitations of the water molecule ||. Especially, one can expect this 
to be the dominant effect in the spectra of multi-bubble sonoluminescence, as the resonances 
discussed under (i) above will average out if bubbles grow and collapse in a random manner. 

5.1.4 Predicted pulse length 

Apart from predicting the radiated spectrum, the theory of quantum vacuum radiation solves 
several conceptual problems that previous theories have not been able to deal with satis- 
factorily. Most importantly, the present theory has no difficulty in explaining the extreme 
shortness of the emitted light pulses; their duration is determined by the parameter 7 de- 
scribing the turn-around time at the collapse and the time it takes for the fluctuations to 
correlate around the bubble. As the latter is on the scale of merely one or a few femtosec- 
onds, 7 is the decisive quantity. Thus one expects the pulse length to lie between lOOfs and 
lOps, which tallies with the experimental observations. 

5.1.5 Absence of radiation in the UV 

Another major question that is successfully answered by the present theory concerns the 
absence of radiation below the absorption edge of water at around 180nm. Water has essen- 
tially no polarizability below this wavelength, so that the real part of its refractive index is 
very close to 1. Hence the mechanism of exciting vacuum fluctuations into real photon pairs 
is inoperative below 180nm; no radiation is emitted and no radiation has to be re-absorbed. 
This explains the absence of any macroscopically discernible effect on the water by the large 
amounts of absorbed light predicted by theories of black-body radiation or bremsstrahlung 
(cf. Sec. I A). 

5.2 Suggested experiments 

Thinking about experiments that distinguish the present from other theories of sonolumi- 
nescence, one quickly comes up with two relatively simple ones. One is to look for pho- 
tons emitted in the X-ray transparency window of water ; both the black-body and the 
bremsstrahlung theories predict a perceptible amount of photons with wavelengths of around 
lA, whereas the present theory denies any photon emission at such short wavelengths since 
the polarizability of water is essentially zero for X-rays, i.e. (n — 1) ~ 0. 

The second presumably easily set up experiment is to force the bubble into an elongated 
rather than spherical shape by using piezoelectric transducers on two or all three axes and 
to examine the angular distribution of the emitted light. For such a case the present theory, 
unlike others, predicts an anisotropic intensity; the number of photons radiated into a given 
direction is roughly proportional to the cross-section of the bubble perpendicular to that 
direction. Thus, if the bubble is spheroidal rather than spherical during the radiation process, 
one expects anisotropy. 
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5.3 Agenda and open questions 

The most important point still to be attended to is to extend the numerical calculations of 
the spectral density (4.5) down to realistic minimum bubble sizes of one micron or less [[IT 



This will allow one to make more precise statements as to the turn-around time 7 required 
to produce the experimentally observed number of photons (cf. Sec. A 1) and to explore the 
effects of resonances if the bubble size is comparable to the photon wavelengths (cf. Sec. A 
3). 

Another effect to be studied in detail is the photon radiation produced by the rapid 
variation of the refractive index of water due to the rapidly varying compression around the 



outside of the bubble ||31|| . Preliminary estimations have shown that this mechanism is of 
secondary importance for the sonoluminescence problem; to understand the principle of it 
might, however, be useful in view of other applications. 

An academic, but nevertheless interesting question to ask is where the photons actually 
are produced. A model calculation for a one-dimensional moving dielectric filling a half-space 
|TT| has indicated that, although the emission of this kind of quantum vacuum radiation 
depends on the existence of a moving interface between two media of different refractive 
indices, the photons do not come directly from the interface but from within a certain vicinity 
of it, as suggested also by physical intuition. However, inconsistencies between where the 
photons are produced and which the support of the radiation pressure on the dielectric is, 
are buried in the assumption of a perfectly rigid dielectric. 

The hydrodynamics of the bubble has been considered as given in the present work; its 
theory has been quite successfully established ||. However, especially in view of relaxing 
restrictions on the velocity and the acceleration of the bubble surface in that theory and 
of scrutinizing the bubble dynamics at the moment of collapse, the important role of the 
back-reaction of the photon radiation onto the bubble should be recognized. The momentum 
loss due to the radiation process might have a significant effect on the bubble. While the 
emission lasts the equations motion of the liquid-gas interface will be supplemented by a 
frictional force which is roughly proportional to the fourth derivative of the velocity of the 
interface, as discussed just below eq. (|4.7|). 

One of the puzzles that remain is why stable single-bubble sonoluminescence is seen 
only in water, although multi-bubble sonoluminescence has been observed in a variety of 
fluids. The present author's conjecture is that the reason for this is buried in the unusual 
properties of gas solubility of water, which conspire with hydrodynamic mechanisms to 
lead to a exceptionally sharp and violent collapse of a driven bubble. In other fluids such 
conditions might be reached at random, but not in a regular fashion to produce radiation 
from stably maintained bubbles. 



5.4 Credo 

To close a conceptual remark might be appropriate. At first sight, the idea that the burst 
of photons seen in sonoluminescence has its origin in the zero-point fluctuations of the 
electromagnetic field might seem utterly strange, as one tends to think of low-energy photons 
emitted from material media as coming from atomic transitions. Pondering on this, one has 
to admit that all we really know is that photons come from some kind of moving charge 
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or, field-theoretically speaking, from the coupling to a fermion field. As we structure our 
thinking, we are most inclined to consider atoms as the basic entities of all materials and try 
to explain all physical phenomena on this basis. However, there is no reason for so doing; 
we are completely at liberty to mentally re-group these charges in a variety of different ways 
and should choose whichever is most appropriate for the problem at hand. In the case of 
sonoluminescence atoms are obviously not the basic entities to be considered, since atomic 
transitions are about a thousand times slower than a sonoluminescence pulse. Here the basic 
structure of the medium with respect to the radiation process is most suitably thought of 
as an assembly of dipoles with a certain dielectric response. This point of view enables 
one to consider the cooperative response of the charges to the zero-point fluctuations of 
the electromagnetic field, and quantum vacuum radiation emerges as a consequence quite 
naturally. 
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APPENDIX A: HAMILTONIAN FOR A DIELECTRIC 
IN UNIFORM MOTION 

The standard way of deriving a Hamiltonian is to proceed from a Lagrangian density. The 
Lagrangian density for a homogeneous dielectric moving rigidly and uniformly is a function 
of the dielectric constant e of the medium and of the velocity (3 of the medium relative to 
the frame of the observer; 

C = C(e,P). 

It is uniquely determined by the following three requirements: 

(i)lia the limit of j3 = it should reduce to the familiar Lagrangian density for a stationary 
dielectric 

£(e,(3 = 0)= 1 -(^-B^ . (Al) 
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(ii) For an optically transparent medium one should recover the Lagrangian density of 
the vacuum, which, by Lorentz invariance, is independent of the velocity (3; 



C(e=l,P) 



4 r ^ r 



-fE 2 



(A2) 



The symbol denotes the field strength tensor of the electromagnetic field, F^ 
W V A^. Its dual is defined F^ = 1/2 e^ a/3 F a/3 . 

(Hi) L must be a Lorentz scalar. The only true scalars that are quadratic in the fields 
and that depend solely on the field strength F^ u and on the four-velocity m m of the medium 
are F^F^, u^ v u a F au , and u^ u u a F av . 

From the above the Lagrangian density is found to be 

1 „ e- 1 

4 



u^u a F av 



(A3) 



Now the Hamiltonian can be derived by going through the canonical formalism 

n = 

5 A 



n 

This leads to the Hamiltonian density 

e 



n-A-c. 



1 6(1 -P") 



n 2 „ 2 

— + B 2 



e-P 
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(A4a) 
(A4b) 

• (A5) 



Substituting II = — D and, with a perturbative treatment in mind, expanding in powers of 
the velocity P, one obtains 



H 



d^r 



+ i_I /5 .( DA B) + 0(/3 2 ) 



(A6) 



This is a very natural result; the Hamiltonian for a stationary dielectric is augmented by an 
energy-fluxlike correction which vanishes for a transparent medium. 

Another way of arriving at the same result is to appeal to the Lorentz invariance of the 
Maxwell theory. The Hamiltonian density must in any frame be given by 



n 



i 



; D E + B H) 



(A7) 



where the fields are as measured in this frame. For a non-magnetic dielectric B = H, but 
the D and E fields are connected by some non-trivial constitutive relation which can be 
found by Lorentz-transforming the constitutive relation D' = eE! valid in the rest-frame of 
the medium into the laboratory frame. There the constitutive relations read 



E|, 
E, 



Dm 



e-P 2 



-(i-P 2 ) r> ± - £ - 



pAB 



(A8a) 
(A8b) 



Utilizing these to replace E and H in eq (|A7|) one recovers the Hamiltonian density ( 
obtained earlier by different means. 
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APPENDIX B: MODE EXPANSION FOR THE HELMHOL' 
EQUATION IN SPHERICAL COORDINATES 



The mode functions AT^ i ^ :m m t ne expansions (|2.7 ) are solutions of the Helmholtz equation 



u 2 e{r) A(r,k) + V 2 A(r, k) = , u = |k| 
Rewriting this equation as 



f^V 2 ^ v^A = -^ 2 v^A 



(Bl) 



(B2) 



makes obvious that this is the eigenvalue equation of the Hermitean operator {1/ y/e) V 2 (l/v / £), 
and hence the mode functions a/^A form a complete set of orthogonal functions. 

In order to diagonalize the Hamiltonian (|2.1| ) into the Hamiltonian ( [2.6|) of the photon 



field by means of the mode expansions (|2.7|), the mode functions should satisfy the orthonor- 
malization conditions 



d r sir] 



A 



l V, (r: k)A™*(r; k') + A™(r; k)A™*(r; k')] = 5 (3) (k - k') (B3) 



(i) 



for each the TE and the TM polarizations. 

To find the solutions of the Helmholtz equation (|B1|) one conveniently proceeds from the 
scalar solution 



Ait 



(2tt) 3 / 2 



$> _w < i e [cosS e j e {nkr) + sin^ y e {nkr)} Y e m *{k)Y e m {r) 



(B4) 



which is normalized to behave like a plane wave e mk ' r /(27r) 3 / 2 for kr — ► oo. The phase Si 
will be chosen to meet the required continuity conditions across the surface of the bubble. 

The vector field operators invariant under rotation are r, V, L = — ir x V, and V x L. 
Since r fails to commute with V 2 , only the last three operators may be used to generate 
rotation-invariant vector solutions of the Helmholtz equation. However, V$ is an irrotational 
field; only the solenoidal fields L$ and (V x L)$ can be employed for representations of the 
transversely polarized electromagnetic field. Choosing Am ~ L$ and A( 2 ) ~ l/(nk) V x L$ 



and observing the correct normalization (|B3[) one obtains for the mode functions outside the 
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bubble 



- E / [cos5 Ul (nA;r)+sin^y,(nA;r)] F £ m *(k) V-/ Y e m (r 

^(i) v = \ ~ E / ; [cos 5 e jt{nkr) + sin <J< ye(nkr)} Y £ m *(k) (-i) 

V 7T £ m , /^Y/> _|_ 1 ^ (7(7 



A (1)r = 



e 


-i5 £ ^ 


V M 


'£+ 1) 


e" 


-iSt j£ 




(* + !) 






e" 




\A< 


'£+!) 


c~ 


-iS e [£ 




(t + 1) 


e" 


1 



/n e~ 1 i 1 

- V , —— [nkr cos 5 t j e {nkr) +nkr sin 8 t y e {nkr)] 

tt 7t JP(P+ -\\ nkr 



dY e m (r) 



x^(k) i 



/n e~ iSe \ l 1 

- V — — [nkr cos 5 e j e (nkr) + nkr sm6 e y ( {nkr)} 
71 tl J 1(1 + 1) nfcr 



Aa)r = \/- E / , t cos ^ £(«**0 + Sin ^ ^ m *(k) Y e m (r) , 

(B5) 



and for those inside the bubble 



^^S S7, 7ra J ' (tr)Fr(k)H) T' 

A(i) r = o , 

^ = ^ E ^ "7=^= T ^ ^ 1 ^W 1 ' (B6) 

V 71 J 1(1 + 1 fcr 99 



>w - f n g «r' -7^^ *w »rw • 



Here and in the following a prime behind a bracket means a derivative with respect to the 
argument of the spherical Bessel function. The mode functions for the inside of the bubble 
have zero phase shift because the fields have to be regular at the origin r = 0, which excludes 
any contributions from the spherical Bessel functions of the second kind ye as these diverge 
for zero argument. 

The phase shifts 5j E '™ and the normalization constants iSj E,TM _1 are determined by 
the continuity conditions fl2.4|) across the bubble surface at r = R. One obtains 

\f TE \f™ 
tan5, TE = |W, tan5™ = ^ f , (B7) 



Sonoluminescence 



24 



and 



where 



<?TE -1 _ _^ * cTM -1 _ _^ * mo\ 

' " (nkR) (-Vj*-iJ^ E ) ' * ~ (kR) (-P™-ijV7 M ) ' 1 J 

A/"™ = ^(jfci?) [nA;/? ^(nibi?)]' - jiinkR) [kR j e (kR)}' , 
VJ E = y e (nkR) [kR j e (kR)}' - 3e (kR) [nkR y e (nkR)}' , 
N™ = 3e (kR) [nkR ji(nkR)}' - n 2 j e (nkR) [kR j e (kR)}' , { > 

V™ = n 2 y e (nkR) [kR Je (kR)}' - Je (kR) [nkR y e (nkR)}' . 



Note that 



lim Sj™ - 1 = 1 , (BIO) 



due to the fact that in this limit the V e ' are simplified by the Wronskian of the spherical 
Bessel functions 

ji(x)y' t (x) -j'e(x)y e (x) = — , 

Ob 



and the J\f^ E '™ obviously reduce to zero. 



APPENDIX C : FORCE ON A STATIONARY DI- 
ELECTRIC 

There are several ways of deriving an expression for the force applied by an electromagnetic 
field on a dielectric body; the physically most intuitive one is to consider the force as ensuing 
from induced currents and surface- charge densities. 

In Lorentz gauge or in Coulomb gauge without free charges the vector potential satisfies 
the wave equation 

| (sA) - V 2 A = . (CI) 

On rewriting this equation as 

A - V 2 A = j ind , 
one obtains an induced current density 

Jind = -(£ - 1)A . (C2) 

By continuity 

dPind v . _ n 
"777 H V ■ Jind — U 

at 

one finds that the induced surface-charge density is 

Pind = -V • ( — D) . (C3) 



£ 

Therefore a stationary dielectric is acted upon by a force density 

f = PindE + jind A B . (C4) 
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On integrating this density over a dielectric with a cavity one has to bear in mind that the 
radial component of the electric field is not continuous across the boundary. However, the 
gradient of (e — l)/e in ( |C3| ) brings about a 5 function on the surface of the cavity, which 
is multiplied by the electric field in eq (|04j) for the force density The mathematically and 
physically correct prescription is to substitute the electric field by the average of its values 
on the two sides of the boundary. Then one obtains for the radial component of the force 
on a spherical bubble of radius R 



i 



i 



n- 



2 



n- 



l + ^)Dl+Bl-Bl 



m 



(C5) 



where the integral runs over the complete solid angle; the tangential components of the force 
are zero as obvious from symmetry. 

Strictly speaking, the expression ( |C5| ) is of course the force on the dielectric and not the 
force on the bubble. Nevertheless, since the two are complementary and the force density 
has 5-support on the boundary, i.e. the force is really applied only to the interface, it seems 
reasonable to speak of the force as acting on the bubble. 

Another, less palpable and more formal way of calculating the force density (|C4|) is 
to proceed from the stress-energy-momentum tensor of the electromagnetic field. In this 
approach care must be taken when interpreting formulae, because the momentum density of 
the photon field in a dielectric medium is subject to ambiguity. However, this issue will not 
be addressed here as the force density is unequivocally defined and interpretable. 

The stress exerted by the fields alone, i.e. the fields as in vacuum and exclusive of the 
polarization fields inside the dielectric, is given by the space-like components of the stress- 
energy-momentum tensor in vacuum 



T^(e = l) = -E l E,-B l B J 



- 5j 
2 ' 



fE 2 + B 



(C6) 



In vacuum there is no doubt about the momentum density carried by the field; it reads 



(C7) 



which is just the Poynting vector. Overall momentum balance requires that the change in 
the mechanical momentum density of the material (i.e. the force density on the dielectric) 
together with the change in the momentum of the fields alone flC7Q equal the negative 
gradient of the stress (|C6|) due to the fields, 



0_ 

dt 



M 



mcch 



-v/^( £ =i; 



Thus the force density is given by 



fi 



dM 



mcch 



Of 



"at 1 ® 



A few trivial transformations taking into account vector identities and the Maxwell equations 
yield 

fi=(l- i) (flfcVfeBi - ~ V.B 2 ' 
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Upon integration over the bubble the first part of this expression is easily seen to lead to 
the -B-dependent terms in the force (|C5|) . In order to recognize the second part one should 
note that VjEj = V • E is zero inside as well as outside the dielectric, but non-zero at the 
interface; Vj-E/j gives a 5 function on the surface multiplied by the difference of the outer 
and the inner electric fields. Thus one recovers the induced surface-charge density ( p3|) , the 
force on which engenders the same _D r -dependent term as in ( |C5l ) before. 
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Figure captions 
Figure 1 

The model profile of eq. (4.9) for the squared radius at the collapse of the bubble, printed 
for four different values of the parameter 7. The solid line corresponds to the function with 
the largest 7. 

Figure 2 

The spectral density calculated numerically from eq. (4.5) as a function of photon frequency, 
for the model profile (4.9) with R = 45/zm, i? m i n = 3/im, and 7 = lOfs. 

Figure 3 

The same data as in figure 2 but as a function of photon wavelength. Note that one can 
barely make out the features that are clearly visible in figure 2. 
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